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1 Summary. Results of 2D simulations of the magnetorotational mechanism of su- 

\Q ' pernova type II are presented. Amplification of toroidal magnetic field of the star 

due to differential rotation of the star leads to the transformation of the rotational 
(gravitational) energy to the energy of the supernova explosion. In our simulation 
the energy of the explosion is 1.12 • 10 51 er<?. The explosion ejects about O.llM©. 
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1 Magnetorotational mechanism 

The magnotorotational supernova (MRS) explosion model was suggested in [1] . The 
idea of MRS consists of getting explosion energy from the rotational (gravitational) 
energy of the collapsed magnetized massive star. ID numerical simulation of the 
C^h' MRS mechanism has been made in [2], [3]. We have made 2D numerical simulation 

J- of the MRS using specially developed implicit conservative Lagrangian scheme on 

triangular grid with grid reconstruction. Our results show that MRS leads to the 
' energy output of the 1.12 • 10 erg and ejection of O.IIMq. 

, Core collapse of a star leads to formation of the rapidly (almost rigidly) rotating 

neutron core and differentially rotating vast envelope. The explosion energy for the 
supernova is taking form the rotational (gravitational) energy of the magnetized 
, star. The magnetic field plays the role of the "transmission belt" for the rotational 

5^ ' (gravitational) energy to the energy of the supernova explosion. Toroidal compo- 

nent of the magnetic field is amplifying with time due to the differential rotation 
of the star. When force produced by magnetic pressure substantially changes the 
local balance of forces a compression MHD wave appears and goes through the 
star's envelope outwards. Moving along steeply decreasing density profile this wave 
transforms to the MHD shock which produces supernova explosion. 
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2 Formulation of the problem 
2.1 Basic equations 

Consider a set of magnetohydrodynamical equations with selfgravitation and with 
infinite conductivity: 

dx dp _ 

du „/ H H\ V-(H<g>H) „^ 
P M=- V ( P+ —) + 4, 

P^-A — )=HVu, Zi<? = 47rGp, 



V P 

p^+ P V -u + pF(p,T) = 0, 

where = + u- V is the total time derivative, x = (r, ip, z), u is velocity vector, p 
is density, p is pressure, H = (H r , H v , H z ) is magnetic field vector, $ is gravitational 
potential, e is internal energy, G is gravitational constant, H ® H is tensor of rank 
2, F(p,T) is the rate of neutrino losses, other notations are standard.. 

Axial symmetry = 0) and symmetry to the equatorial plane [z = 0) are 
assumed. 



2.2 Equations of state 

The equations of state P(p,T) [4]. includes approximation of the tables from [6, 8] 
for the cold degenerate matter P (p): P = P(p, T) = P (p) + pKT + 2^1, 

PM = / P o (1) =&ip 1/3 /(l + cip 1/3 ), a*p< P i, 

bi = 10.1240483, a = 10~ 2 ' 257 , pi = 10 9 ' 419 , b 2 = 1.0, c 2 = 1.1598, p 2 = 
10 11 ' 5519 , b 3 = 2.5032, c 3 = 0.356293, p 3 = 10 12 ' 26939 , 6 4 = 0.70 40 15 15, c 4 = 
2.117802, p 4 = 10 14 ' 302 , 6 5 = 0.16445926, c 5 = 1.237985, p 5 = 10 15 0388 , & 6 = 
0.86746415, c 6 = 1.237985, p 6 > p 5 , a = 10 26 1673 ' 

Here p is a total mass-energy. The energy of the unit mass is defined as: e = 
£o(p) + §KT + , where 5ft - gas constant, a - constant radiation density, and 
p 



eo( P ) = J^dp. 


The neutrino losses from Urea processes are defined by relation [7]: 

_ 1.3 ■ 10 9 ae(T)r 6 _ x _j 

f(P,T) = ^— r erg-g -c , 2) 

1 + (7.1 ■ io- 5 pr) ! 

1, T< 7, 

ae(T) = { 664.31 + 51.024(T - 20), 7 < T < 20, (3) 

664.31, T > 20, T = T ■ 10~ 9 . 
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Neutrino losses from pair annihilation, photo-production and plasma were also taken 
into account. These type of the neutrino losses have been approximated by the 
interpolation formulae from [9] : Q to t = Q V a%r + Qphoto + Qpiasm These three terms 
can be written in the following general form: 

Qd = ^' Q)e eTMTMW' (4) 

For d = pair, K(p,a) = g(a) e -^, g(a) = 1 - *|§4 + ^ + ^ + For d = 

photo, K(p,a) = (p/nz)a~ 5 ; For d = plasm, K(p,a) = (p/pz) 3 ; £ = { E {§f-) a - 
Here pz = 2 is number of nucleons per electron. Coefficients c, ai, b» for different d are 
given in [9]. The general formula for neutrino losses in nontransparent star has been 
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Fig. 1. Initial poloidal magnetic field. 



written in the following form: F(p,T) — (f(p,T) + Q to t)e~T^. The multiplier e~TO , 
where t v = S v nl v restricts neutrino flux for non zero depth to neutrino interaction 
with matter r„ . The cross-section for this interaction S v was represented in the form 
S v = ,„ Trrnm , the nucleons concentration is: n = — , m v = 1.67 • 10 _24 g. The 
characteristic length scale l v which defines the depth for the neutrino absorbtion, 
was taken to be equal to the characteristic length of the density variation, as l v = 
T=£-r = TT7Q-- This value monotonically decreases when moving to 

\ V P\ ((a P /Sr)2 + (9p/Sz) 2 ) 1 

the outward boundary, its maximum is in the center. It approximately determines 
to the depth of the neutrino absorbing matter. The multiplier 1/10 was applied 
because in the degenerate matter of the hot neutron star only part of the nucleons 
with the energy near Fermi energy (it was taken w 1/10) takes part in the neutrino 
interaction processes. 



3 Results 

For the numerical simulations we have used implicit Lagrangian difference scheme on 
triangular grid with grid reconstruction. For the description of the applied numerical 
method see, for example, [5] and references therein. The number of knots of the 
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triangular grid was about 5000. As a first stage of the MRS mechanism we have 
calculated a collapse of rotating star, leading to differentially rotating configuration. 
After the collapse the star (core of evolved massive star) consists of a almost rigidly 
rotating neutron star with radius ~ fOkm which rotates with the period ~ O.OOf 
sec, and vast differentially rotating envelope . After formation of the differentially 




Time (sec.) Time (sec.) 

Fig. 3. Time evolution of the (ejected energy) /(rotational energy) in % (left plot) 
and (ejected mass)/(total mass) in %. (right plot). 

rotationg configuration the initial poloidal magnetic field was switched on (Fig.f). 
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The energy of the initial magnetic field was E ma gO = W~ 6 E gr av Where E mag o 
- is the energy of the initial magnetic field, E gra v - is the gravitational energy of the 
collapsed star. 

Due to the differential rotation the toroidal component of the magnetic field 
is increasing with time. The magnetic pressure growths and produces compression 
MHD wave, which moves through envelope with steeply decreasing density. Soon 
after appearing this wave transforms into the MHD shock, which throws away part 
of the matter of the envelope. The MHD shock front is clearly seen at the velocity 
field plot (Fig. 2 - left plot). The magnetic field transmits angular momentum of the 
neutron star outwards (Fig. 2 - right plot). 

Results of our simulations show that the energy of the supernova explosion is 
about 1.12 • I0 51 erg (35% of the rotational energy of the star). The explosion ejects 
about O.IIMq (~ 9.7% of the mass of the star). At the Fig. 3 the time evolution of 
the ejected mass and ejected energy are presented. 

Detailed description of the results of the simulations of the MRS will be published 
elsewhere. 

The authors G.S.B.-K. and S.G.M would like to thank the Organizing Committee 
for the support and hospitality. This work was partially supported by the grant 
RFBR 02-02-16900. 
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